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The isoscalar giant monopole resonances (ISGMR) and giant dipole resonances (ISGDR) in medium-heavy nuclei 
are investigated in the framework of HF+RPA and HF-BCS+QRPA with Skyrme effective interactions. It is 
found that pairing has little effect on these modes. It is also found that the coupling of the RPA states to 2p-2h 
configurations results in about (or less than) 1 MeV shifts of the resonance energies and at the same time gives the 
correct total widths. For the ISGMR, comparison with recent data leads to a value of nuclear matter compression 
modulus close to 215 MeV. However, a discrepancy between calculated and measured energies of the ISGDR in 
^"**Pb is found and remains an open problem. 



Introduction 

The determination of the nuclear incompressibihty K is 
still a matter of debate, despite a remarkable number 
of works on the subject'^-'. In the present contribution, 
we present self-consistent calculations of the nuclear col- 
lective modes associated with a compression and expan- 
sion of the nuclear volume, namely the isoscalar giant 
monopole and dipole resonances (ISGMR and ISGDR, 
respectively). In fact, we share the point of view that 
the most reliable way to extract information on K is to 
perform that kind of calculations, having as the only phe- 
nomenological input a given effective nucleon-nucleon in- 
teraction, and choose the value of K corresponding to the 
force which can reproduce the experimental properties of 
the compression modes in finite nuclei. 

The ISGMR, or "breathing mode" , is excited by the op- 



erator 



and it has been identified in many iso- 



topes along the chart of nuclei already two decades ago. 
However, this systematics has never allowed an unam- 
bigous determination of This was one of the moti- 
vations for the recent experimental program undertaken 
at the Texas A& M Cyclotron Institute, which has al- 
lowed the extraction of experimental data for the ISGMR 
of better quality as compared to the past, by means 
of the analysis of the results of inelastic scattering of 
240 MeV a-particles. We refer to other contributions 
in these proceedings for reports on these experimental 
data^-* . Monopole strength functions turn out to be quite 
fragmented for nuclei lighter than ^"Zr. For nuclei like 
208pb, "^Sm, "6sn and ^"Zr, however, one is able to 
identify a single peak which, together with a high-energy 
extended tail, exhausts essentially all the monopole En- 
ergy Weighted Sum Rule (EWSR). These medium- heavy 
nuclei are, therefore, those suited for the extraction of 
information about the nuclear incompressibihty and we 
concentrate ourselves on them in the present work. 



The ISGDR is excited by the operator rf Fio and cor- 
responds to a compression of the nucleus along a defi- 
nite direction, so that it has been called sometimes the 
"squeezing mode" . Although some first indication about 
the energy location of this resonance dates back to the 
beginning of the eighties, a more clear indication about 
its strength distribution in ^osp]^ j^g^g j^gen reported only 
recently"^^. Measurements have been done also for other 
nuclei, namely ®°Zr, ^^^Sn and ^''^Sm (as in the case of 
the giant monopole resonance). There is some expecta- 
tion that the study of this mode can help to shed some 
light on the problem of nuclear incompressibihty. Ac- 
tually, at first sight this compressional mode seems to 
provide us with a new problem. A simple assumption 
like the scaling model (illustrated for the present pur- 
poses in Ref.*^) would lead to two different values of the 
finite-nucleus incompressibihty Ka if applied to the IS- 
GMR and the ISGDR with the input of their experimen- 
tal energies. The hydrodynamical model gives two results 
which are closer'^^ but which still make us wonder about 
the validity of methods based on extracting Ka and ex- 
trapolating it to large values of A, for the determination 
of K. This points again to the necessity of reliable micro- 
scopic calculations of the compressional modes, in order 
to reproduce the experimental data and extract the value 
of K from the properties of the force which is used. 

Our calculations are performed within the framework 
of self-consistent Hartree-Fock (HF) plus Random-Phase 
Approximation (RPA) . We use effective forces of Skyrme 
type'^"''^ and we look at their predictions for the prop- 
erties of ISGMR and ISGDR. The parametrizations we 
employ span a large range of values for K (from 200 MeV 
to about 350 MeV). In particular, we focus mainly on two 
original aspects: firstly, we look at the effects of pair- 
ing correlations in open-shell nuclei; secondly, we study 
if the picture obtained at mean-field level is altered by 
the inclusion of the coupling of the giant resonances to 
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more complicated nuclear configurations. This inclusion 

is necessary if one wishes to understand theoretically all 
the contributions to the resonance width and may shift 
the resonance centroid. 

About the first aspect, it is well known that pairing corre- 
lations are important in general to explain the properties 
of ground states and low- lying excited states in open-shell 
nuclei. Since we wish to see how these correlations affect 
in particular the compressional modes, we take them into 
account by extending the HF-RPA approach to a quasi- 
particle RPA (QRPA) on top of a HF-BCS calculation. 

About the second aspect, we recall that if we start from 
a description of the giant resonance as a superposition of 
one particle-one hole (Ip-lh) excitations, in their damp- 
ing process we must take care of the coupling with states 
of 2p-2h character. They are in fact known to play a 
major role and give rise to the spreading width of the 
giant resonance which is usually a quite large fraction 
of the total width. Within mean field theories, only the 
width associated with the resonance fragmentation (Lan- 
dau width) and the escape width are included (the lat- 
ter, provided that Ip-lh configurations with the particle 
in the continuum are considered). In the past, we have 
developed a theory in which all the contributions to the 
total width of giant resonances are consistently treated 
and we have obtained satisfactory results when applying 
it to a number of erases. In particular, we will recall what 
has been obtained*^' for the case of the ISGMR in ^"Spb. 
We also report about a new calculation for the ISGDR 
in the same nucleus. 



Formalism: a brief survey 

For all nuclei we consider, we solve the HF equations on 
a radial mesh and, in the case of the open-shell isotopes, 
we solve HF-BCS equations. A constant pairing gap A 
is introduced (for neutrons in the case of ^^^Sn and for 
protons in the case of ^°Zr and ^''^Sm), and at each HF 
iteration the quasi-particle energies, the occupation fac- 
tors and the densities to be input at the next iteration are 
determined accordingly. A is obtained from the binding 
energies of the neighboring nuclei^^ . The states included 
in the solution of the HF-BCS equations are those be- 
low a cutoff energy given by Xhf + 8.3 MeV {Xhf being 
the HF Fermi energy), in analogy with the procedure of 
Ref.io). 

Using the above self-consistent mean fields we work out 
the RPA or QRPA equations (respectively on top of HF 
or HF-BCS), in their matrix form. Discrete positive en- 
ergy states are obtained by diagonalizing the mean field 
on a harmonic oscillator basis and they are used to build 
the Ip-lh (or 2 quasi-particles) basis coupled to J'^=0~^ 
or 1~ . The dimension of this basis is chosen in such a way 
that more than 95% (typically 97-99%) of the appropri- 



ate EWSR is exhausted in the RPA or QRPA calculation. 

More details, especially on the way the QRPA equations 
are implemented, will be given in Ref.^^-*. 

As mentioned in the previous section, in the case of ^°*Pb 
we perform also calculations that go beyond this simple 
discrete RPA. This is done along the formalism described 
in Ref.^^^, which is recalled here only very briefly. 

We label by Qi the space of discrete Ip-lh configurations 

in which the RPA equations are solved. To account for 
the escape width F^ and spreading width F-^ of the giant 
resonances, we build two other orthogonal subspaces P 
and Q2- The space P is made of particle-hole configura- 
tions where the particle is in an unbound state orthogonal 
to all the discrete single-particle levels; the space Q2 is 
built with the configurations which are known to play a 
major role in the damping process of giant resonances: 
these configurations are Ip-lh states coupled to a collec- 
tive vibration. Using the projection operator formalism 
one can easily find that the effects of coupling the sub- 
spaces P and Q2 to Qi are described by the following 
effective Hamiltonian acting in the Qi space: 



n{E) = QiHQi + W^iE) + W^{E) 

= Q,HQ, + Q,HP^-^L__.PHQ, 

E - Q2HQ2 + le 

where E is the excitation energy. For each value of E 
the RPA equations corresponding to this effective, com- 
plex Hamiltonian 'H{E) are solved and the resulting sets 
of eigenstates enable us to calculate all relevant quanti- 
ties, in particular the strength function associated with a 
given operator. To evaluate the matrix elements of W^, 
we calculate the collective phonons with the same effec- 
tive interaction used for the giant resonance we are study- 
ing (within RPA), and we couple these phonons with the 
Ip-lh components of the giant resonance by using their 
energies and transition densities. 



Results for the isoscalar monopole resonance 

As recalled in the introduction, Youngblood et al?^ have 
recently measured the ISGMR strength distribution with 
fairly good precision, in the nuclei ^°Zr, ^^^Sn, ^^^Sm and 
^"^^Pb. In their work, they also compare the experimental 
centroid energies with the calculations of Blaizot et al}^'> 
performed by using RPA and employing the finite-range 
Gogny effective interaction: a value of the nuclear incom- 
pressibility K = 231 MeV is deduced. In the following, we 
denote as centroid energy the ratio Eq = mi /mo (mo and 
mi being the non-energy-weighted and energy-weighted 
sum rules, respectively). 

If we try to compare the experimental values with calcu- 
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lations done at the same RPA level but using the zero- 
range Skyrme effective interactions, wc can infer a differ- 
ent conclusion with respect to the value of K. Among 
the Skyrme type interactions, the parametrization which 
gives probably the best account of the experimental cen- 
troid energies in the nuclei studied by the authors of 
Ref.2), is the SGII force^). The results are shown in Ta- 
ble 1. The force SGII is characterized by a value of the 
nuclear incompressibility K = 215 MeV: since it repro- 
duces very well the ISGMR centroid energy in ^°^Pb, and 
it slightly overestimates those in the other isotopes, one 
would conclude that K is of the order of or slightly less 
than 215 MeV. 

This conclusion is inferred by means of simple RPA. It is 
of course legitimate to wonder if calculations beyond this 
simple approximation could lead to different values of the 
nuclear incompressibility. We first consider the effect of 
pairing correlations. In the case of ^^^Sn, the centroid 
energy of 17.18 MeV obtained with the SGII force in 
RPA, becomes 17.19 MeV if one turns to QRPA. A very 
small shift is found also when other forces are used (for 
instance, with the recently proposed SLy4 forcc^', one 
obtains 17.51 MeV and 17.59 MeV for RPA and QRPA, 
respectively) and when other nuclei are considered. In 
general, although wc know that pairing correlations play 
a crucial role not only to explain the ground-state of 
open-shell nuclei but also their low-lying excited states, 
it appears that they do not affect so much the giant res- 
onances like the ISGMR (or ISGDR, anticipating results 
of the next section) which lie at relatively high excita- 
tion energy compared to the pairing gap A. Civitarese 
et al}^^ found also small shifts (of the order of 100-150 
keV) for the ISGMR and ISGQR when pairing correla- 
tions are taken into account: this shift is larger than that 
obtained in the present work, but it is the result of a dif- 
ferent (non self-consistent) model. 

The present conclusion for the nuclear incompressibility 
is therefore similar to that obtained by Hamamoto et 
al.^^\ since they find that the Skyrme interaction which 
provides the best results for the ISGMR is the SKM* 
parametrization and this is very similar to SGII (the as- 
sociated nuclear incompressibility being 217 MeV). Our 
study, however, is done in a more general framework since 
we have analyzed also the role of pairing correlations. 

If we finally consider the results of calculations beyond 
mean ficld^-' (which include not only the continuum cou- 
pling but also the coupling with the 2p-2h type states) 
performed for ^°^Pb we find that it is also possible to re- 
produce rather well the total width of the ISGMR, which 
is around 3 MeV. This width is actually in large part a 
consequence of fragmentation (or Landau damping): at 
least three states share, at the level of RPA, the resonance 
strength, but continuum as well as 2p-2h couplings are 
able to give to each peak the correct width so that the 
overall lineshape coincides with the experimental find- 
ings. We stress that the coupling with the 2p-2h type 



states is also responsible for a downward shift of the IS- 
GMR centroid and peak energies, which is of the order of 
0.5 MeV. One may argue that this affects the extraction 
of the value of K from theoretical calculations. Actu- 
ally, since the value of K associated with a given force is 
obtained by a calculation of nuclear matter at the mean 
field level, it is legitimate to draw conclusions about K 
from the comparison with the experiment of the ISGMR 
results for finite nuclei obtained again at the mean field 
level. But the fact that a given force is able to account for 
the ISGMR linewidth enforces our confidence about its 
reliability. And it would be of course legitimate either, to 
compare the centroid energies obtained after 2p-2h cou- 
pling with experiment provided the value of K associated 
with the force is calculated by including the same cou- 
plings at the nuclear matter level. No such calculations in 
nuclear matter have been done so far, to our knowledge. 

Table 1. Experimental and theoretical values of the cen- 
troid energies Eq = mi /mo for the ISGMR and IS- 
GDR. The theoretical values are obtained with the 
Skyrme-RPA approach, using the SGII interaction. 
All values are in MeV. 





ISGMR 


ISGDR 


Exp. 


Theory 


Exp. 


Theory 




17.9 


19.1 


26.2 


27.1 




16.0 


17.2 


23.0 


26.3 


i44Sni 


15.3 


16.2 


24.2 


25.4 


208pb 


14.2 


14.1 


20.3 


24.1 



Results for the isoscalar dipole resonance 

A peculiar feature of calculations of this giant resonance 
is the appearance of a spurious state in the calculated 
spectrum. When diagonalizing the RPA matrix on a 
1~ basis, we expect to see among all states the spuri- 
ous state at zero energy corresponding to the center-of- 
mass motion and we expect as well that it exhausts the 
whole strength associated with the operator X^i'^i^io- 
Due to a lack of complete self-consistency (some part 
of the residual interaction, like the two-body spin-orbit 
and Coulomb forces, are usually neglected in the RPA 
because their effect should be rather small) and to nu- 
merical inaccuracies, this is not the case. The spurious 
state comes out in practice at finite energy and its wave 
function does not overlap completely with that of the 
exact center-of-mass motion: as a consequence, the re- 
maining RPA eigenstates are not exactly orthogonal to 
the true spurious state, and their spurious component 
must be projected out. This is not difficult if RPA is 
done in the discrete p-h space. In any case, it can be 
shown^^^ that this projection procedure is equivalent to 
replacing the J2- rf Yio operator with ^^irf — ?7r,:)Yio, ij 
being |(r^). Once this projection is done, we find that a 
substantial amount of r^Yio strength still remains in the 
10 - 15 MeV region, in addition to the strength in the 
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20 - 25 MeV region. The data do not show any strength 

in the lower region, i.e., experimental centroid energies 
correspond to the energy region above 15 MeV. There- 
fore, to have a meaningful comparison with experiment 
we will refer from now on to theoretical centroid energies 
calculated in the interval 15 - 40 MeV. 

In Table 1 we show the ISGDR centroid energies obtained 
with the SGII force. Especially in the case of ^°®Pb and 
^^^Sn, it can be noticed that RPA calculations tend to 
overestimate the value of the centroid energy, the dis- 
crepancy being less severe in the other two cases. One 
may wonder if this is a special feature of SGII, although 
this force has been said to behave rather well for the 
monopole case. Fig. 1 shows that this is not the case: the 
centroid energies obtained in RPA with a number of dif- 
ferent Skyrme parametrizations are plotted as a function 
of their incompressibility if, and it can be noticed that 
all forces systematically overestimate the experimental 
values of the centroid energies. Gogny interactions have 
also been used to study the ISGDR in this and other 
nuclei^ ''^ but they also predict too large centroid ener- 
gies. The same can be said about relativistic models like 
relativistic RPA^^^ or time-dependent relativistic mean 
fieldi9). We may conclude that the case of the ISGDR in 
208p^ is a kind of exception among the giant resonances 
studied within the self-consistent HF-RPA approach, as 
usually one never finds such large discrepancies between 
theory and experiment (~ 4-5 MeV). 

We finally address the question whether the coupling 
with more complicated configurations, which has been 
seen to be responsible of a downward shift of the reso- 
nance energy, can diminish this discrepancy in the case of 
208p|-, -y^g have done for the ISGDR a calculation of the 
type described above in the monopole case. The resulting 
strength function, which includes continuum and 2p-2h 
coupling, is depicted in Fig. 2. One can see that the 
total width of the resonance is quite large, and the theo- 
retical value of about 6 MeV compares well with the ex- 
perimental result which is about 7 MeV^°) . Although the 
resonance lineshape is accounted for by theory, the down- 
ward shift with respect to the RPA result is only about 1 
MeV. The fundamental problem why the ISGDR energy 
in ^°^Pb cannot be reproduced by theoretical models still 
remains. 



Conclusion 

In this paper, we have considered the isoscalar monopole 
and dipole resonances in a number of nuclei and we have 
tried to reproduce their properties by means of HF-RPA, 
or QRPA on top of HF-BCS, or more sophisticated ap- 
proach which takes care of the continuum properly and of 
the coupling with nuclear configurations which are more 
complicated than the simple Ip-lh. In general, we have 
found that the effect of pairing correlations is quite small 



as these resonances lie at high energy with respect to the 
pairing gap A. 

Concerning the RPA results, the situation looks differ- 
ent in the case of monopole and dipole. In the former 
case, the Skyrme-type force SGII is able to reproduce 
well the centroid energy in ^°^Pb, and it slightly overes- 
timates this energy in other medium-heavy nuclei which 
have been measured accurately in recent experiments. 
This would allow us to extract a value of the nuclear 
incompressibility around 215 MeV. In the case of the 
ISGDR, however, the same Skyrme force overpredicts 
this centroid energy in ^°^Pb by about 4 MeV. Other 
parametrizations of Skyrme type cannot do better, and 
the problem is not solved if one turns to Gogny inter- 
actions or to relativistic models. Therefore, although in 
other nuclei this discrepancy between theory and experi- 
ment can be less than in the case of ^°*Pb, we can say that 
the "squeezing mode" , which could be taken as a further 
probe of the nuclear incompressibility besides the well- 
known "breathing" monopole oscillation, is challenging 
us with a new problem. 

Calculations beyond mean field do not change substan- 
tially our conclusions about the centroid energies. How- 
ever, we stress that these calculations are necessary for a 
proper account of the giant resonances lineshape, and in 
fact in our case they have been able to reproduce the to- 
tal width of the ISGMR in ^"^Ph, and also of the ISGDR 
although its centroid is overestimated. 

We would like to thank Uniesli Garg for stimulating dis- 
cussions and for communicating experimental data prior 
to publication, and also Jean-Paul Blaizot for useful dis- 
cussions. 

Fig. 1. RPA centroid energies of ISGDR calculated in 
208 with various Skyrme interactions, as a function 
of compression modulus K. 

Fig. 2. Calculated ISGDR strength distribution of oper- 
ator (r^ — r]r)Yio including damping effects of parti- 
cle-hole-plus-phonon coupling (see text). The nucleus 
is 208pb. 
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